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ABSTRACT 


PUMPING-EROSION SUBSIDENCE 
OF A SEAFLOOR PLATE-FOOTING 


(November 1987) 


ALVIN EUGENE GRIMMIG JR., B.S. University of Washington 


Chairman of Advisory Committee: Dr. John B. Herbich 


Subsidence of objects that rest on the seafloor is an 
expected but not well understood phenomenon. Subsidence is 
the settlement of fixed structures and mobile structures 
(such as oceanographic packages) that are placed on the 
seafloor. A major concern involving bottom resting systems, 
in relatively shallow water (less than 200 meters) is wave 
induced motions. The system response of the structure in 
this dynamic environment can create substantial loadings at 
the structure-seafloor interface. Past works have addressed 
wave-induced and current scour as the major factors causing 
structural subsidence, neglecting the effects of pumping - 
erosion. 

The research described is directed towards indentifying 
and understanding the phenomenological aspects of pumping- 
erosion. The work describes a new hypothesis concerning 
liquefaction as being the precursor to subsidence in a two 


cycle pumping-erosion process. Due to the complexity of 
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structure-soil interaction a multivariable linear regression 
was performed on the data to substantiate experimental 


observations. 
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CHAPTER I 
INTRODUCTION 


Structural subsidence can range from a few inches to 
several feet causing results varying from slight tilting and 
displacements of a structure to total and catastrophic 
structural failure. As a result of the petroleum and gas 
industries use of seafloor supported structures (i.e. jack- 
up rigs) much investigation has been directed at the 
structure-soil system and interaction. Although this area 
has been intensely investigated, the focus of attention has 


been on large oil rig structures. 


The use of relatively small mobile submerged structures 
for oceanographic and industrial purposes has increased in 
the past several years. The ability to successfully deploy, 
Operate and retrieve these instruments and equipment 
packages is largely dependent upon the structure-seafloor 
interaction. Bottom resting systems are subject to the 
dynamic response of the structure in the ocean regime which 
Can produce significant dynamic loadings at the seafloor. 
The stability of a submerged seafloor structure is dependent 
upon the subsidence and orientation of the structure while 
subject to wave and current forces. These forces produce a 


footing response in a variety of bearing, tension, sliding, 
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rocking and torsional strains in the marine soil supporting 


the structure (Figure 1). 


The central issue of this report involves the 
investigation of the phenomenon and mechanisms of pumping- 
erosion and its contribution to subsidence. Pumping-erosion 
is inherently non-linear as the structure is intermittently 
coupled and decoupled from the soil interface. In addition, 
the effects of the period of oscillation, amplitude of 
displacement, footing characteristics, hydrodynamic effects 
and soil parameters are simultaneously involved resulting in 


a very complicated and complex problem. 


The basic pumping-erosion cycle involves an uplifting 
force vertically displacing the footing off the seafloor. 
Subsequently, a downward force will dynamically load the 
structure onto the soil again. This over simplified cycle 
is a Trojan horse, as it involves several discrete and 


intertwined forces which are not readily apparent. 


Consideration must be given to the geometry and 
dynamics of the footing, The size, shape, and weight of the 
structure will determine such characteristics as center of 
mass, center of buoyancy, drag and lift forces and moments 
Of inertia which will determine the structures response to 
the exciting force. The wave and current environment will 


produce the exciting force characteristics of period, 
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amplitude, and duration of excitation. Hydrodynamically, 
effects of exit and entrance velocities, water density and 
viscosity as well as the added mass of the interacting 
structure will influence the structure response and the 


sediment transport and settlement. 


Soil parameters such as soil density, permeability, 
void ratio, grain size and shape, and pore water pressures 
will determine seafloor response to the structure-soil 


interaction. 


Traditionally subsidence was attributed to the effects 
of scour resulting from the driving force of wave and 
currents with little regard to pumping-erosion. Any 
inference to pumping-erosion usually surrounds the latter 


half cycle where the water is displaced or "pumped" from 
under the footing as it makes contact with the soil surface. 
The soil transport has been attributed to the "jetting" 


action of the water escaping from under the footing. 


But what of the first half of the cycle of the footing 
when it is uplifted from the seafloor ? What forces are 
involved ? And what contributions do these forces make to 
the process of subsidence ? As the footing is uplifted, 
negative pore water pressures can be developed in the soil 


under the footing, possibly to the point of liquefaction. 
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Liquefaction, or fluidization, occurs when the pore- 
water pressure in a soil is equivalent to the total stress 
of the soil; in other words when the effective strength of 
the soil equals zero. Liquefaction of a soil requires 
considerable pore water pressures as compared to the 
relative ease of fluidization, which reduces the effective 
strength of the soil by the increase in the pore water 
pressure combined with the uplifting effect of the upward 
seepage flow. If the bed is considered deformable, the 
fluidization would then weaken the grain skeleton and 
distort the grain bed surface causing the grains to have a 
greater hydraulic profile making them readily available for 


transport during the downward half cycle of the footing. 


Thus, it is hypothesized that subsidence is a two phase 
cycle. The first phase occurs on the uplifting of the 
footing causing sufficient negative pore water pressures to 
fluidize the bed which weakens and distorts the bed making 
the upper grain layers more susceptible to sediment 
transport. The second phase, occurs on the downward stroke 
of the footing where the velocities developed under the 
footing transport the sediment. The transport to the 
sediment will be accentuated at the edges where the water is 
jetted out from under the footing. The following research 
is designed to investigate the hypothesis on a 


phenomenological basis. 
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OBJECTIVES 


The objectives of this research are as follows: 


1. To determine which geometric, 
hydrodynamic and kinematic parameters are involved 


in pumping-erosion subsidence. 


2. 10 evaluate the relative importance of 
the geometric, fluid and kinematic parameters 


involved in pumping-erosion subsidence. 


3. To determine if sufficient negative pore 
water pressures can be developed under the footing 
to fluidize the foundation bed to the point of 


liquefaction. 
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CHAPTER II 
LITERATURE REVIEW 


There are a limited number of references dealing with 
soil/water/structure interaction. Even fewer address the 
soil/water/structure interactions and forces involved in 
structural subsidence. Furthermore, there are only two 
sited works that discuss the topic of pumping-erosion. Most 
structural subsidence is attributed to he effects of scour. 
Herbich et al.(3) has compiled an extensive biography on the 


topic of conventional scour. 


Reimnitz and Kempema (10) in 1981, describes the 
effects of dynamic ice wallowing in the Alaskan nearshore 
areas. The report hypothesizes the occurence oftwo types of 
hydraulic processes. In the first process the ice plays a 
passive role, acting as a flow obstacle in the current and 
wave regime. This particular aspect would account for the 
more conventional scour around structures. The second 
hydraulic process involves ice playing an active role, 
either by simple vertical oscillations or by wallowing 
(rocking). Dahlberg (1) discusses pumping erosion as it 
pertains to structures. He states that pumping erosion is 
associated with excess pore-water pressures set-up in the 


foundation soil during storm periods. Dalhberg suggests 
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that excess pore-water pressures develop to balance the 
overturning moments on the base of the structure. It is 
further stated that if there is free communication between 
the foundation soil and the seabed, the pore-water gradient 


may be high enough to liquefy the soil locally. 


It was verified that a North Sea platform, Frigg CDP1, 
which had no protective scour skirts, was subject to 
pumping-erosion. Divers observed periodic puffs of sediment 
around the periphery of the structure base when the 
sediments were carried in suspension as the structure rocked 
in the storm environment. It is noted that the pumping- 
erosion described was solely attributed to the positive 
pore-water pressures escaping above the mudline due to the 


pressure gradient. 


An extreme case of pumping-erosion involved the 
Christchurch Bay Tower located in 8.4 meters of water. The 
platform was subject to severe storm conditions. Due to the 
absence of protective scour skirts the platform was 
undermined by scour, which lead to a free rocking motion of 
the structure, resulting in pumping-erosion on the 


foundation soil. 


Pumping-erosion of gravity structures is now rare,due 


to their size and the use of scour skirts. Most settlement 
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of gravity structures, as described by DiBiagio (2), is 


attributed to consolidation settlement. 


Teramoto et al. (13) discussed the scour encountered with 
various sit-on-bottom type of structures. The study 
investigates the various scour patterns and contours for 


various footing designs and configurations. 


However, pumping erosion is more applicable to 
lightweight mobile structures (i.e. oceanographic research 
instrumentation packages which can be deployed from research 


vessels) and is the focus of this topic. 


Several works have dealt with wave induced pore water 
pressures. Approaches have varied from considering the soil 
skeleton to be rigid, to being compressible and pore water 
as being compressible or incompressible. Governing 
equations were usually derived from Darcy's Law of flow 
through a permeable bed or Biot's Law for three dimensional 
consolidation for a poro-elastic material. In addition, 
various studies have evaluated whether residual pore-water 
pressures exist in cohesionless soils or whether the pore 


pressures are merely transient. 


Oldinziel and Brink (9), indicate that an upward flow 
of water through a porous sand bed reduces the apparent 


weight of the sand particles and therefore reduces the sand 
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particles' stability. The results of the study conclude 
that upward flow (blowing) through the bed increases the 
rate of sediment transport. The study relates the upward 
blowing forces to pressure gradients in the bed that depend 


on seepage velocities as defined by Darcy's Law. 


Martin and Aral(5), 1971, indicates that the seepage 
force on a surface grain is only 50% of that for an embedded 
grain. However, the study concludes that it is clear that 
upward seepage will reduce the stability of surface grains, 
whereas the downward seepage will increase stability. 

Martin (4) in 1974, demonstrated that bed failure could also 


result from a horizontal pressure gradient. 
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CHAPTER III 
PROCEDURES AND INSTRUMENTATION 


When non-fixed seafloor systems are subjected to wave 
and current forces of sufficient levels, the structure will 
oscillate. The wave pressure forces, in conjunction with 
the current forces will lift the structure and then will 
load it again on the seafloor. The response of the 
structure and subsequent subsidence are a function of 


several geometric, hydrodynamic and kinematic factors. 


Due to the amount of limited data concerning pumping- 
erosion subsidence a simple plate type footing employing 
purely vertical oscillations was chosen. The tests were 
conducted in a steel tank with a glass observation window. 
The dimensions of the tank were 8 feet long, 4 feet wide and 
4 feet deep. The glass observation window was 5/7 inches 
‘long and 3/4 inches thick. A sand bed, 11.5 inches deep was 
iphaced in the steel tank. The test section (4 feet by 4 
feet) was 11.5 inches deep. The remainder of the tank had 


5.5 inches of sand overlying a false bottom (Figure 2). 


The sand utilized for the test was an Ottawa sand with 


a specific gravity of 2.65 and a mean grain diameter (d50) 


of .22mm (0.087 inches) (Figure 3). The sand had a static 


(pluivated) submerged angle of repose of 16.50. 
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Two steel plates were utilized in the experiments. A 
12 inch square plate (1 square foot), 3/4 inch, 25 pound 
plate was utilized for experiments 1 through 6. A 10 gauge 
(0.132 inch) thick, 16.97 inch square (2 square feet) plate 
attached to the 1 square foot plate was utilized for the 


experiments 7 through 12. 


The plate was uplifted by a cable and pulley system 
coupled to a 1/2 horsepower D.C. variable speed motor via a 
rotating cam follower that engaged lever arm attached to the 
cable (Figure 4). This particular arrangement allowed the 
plate to be uplifted from the sandbed to a predetermined 
amplitude and then the plate was released and allowed to 
free fall onto the test bed. The next cycle would engage 
the cam follower on the lever arm lifting the plate again, 
repeating the cycle. The plate was thus able to follow the 


contours of the scour caused by the pumping-erosion. 


The amplitude of oscillation (1 or 0.8 inches) was 
obtained by eccentrically locating the cam follower from the 
center of rotation of the motor shaft. The selection of the 


amplitudes was restricted by equipment limitations. 


Oscillation periods of 5, 10 and 20 seconds were 


selected to simulate general storm wave conditions that 


might be experienced by seafloor systems. 
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Pore-water pressures in the test bed were measured by a 
Statham P-22S-350 bi-directional pressure transducer (2 
psid). The transducer was placed outside of the tank at the 
same equivalent hydrostatic depth as the sand bed. The 
transducer detected pressure fluctuations through the active 
port via a 1/8 inch clear PVC tubing coupled to a needle 
shaped copper probe fitted with a # 200 mesh screen to 
restrict the entrance of sand particles into the transducer 
port. The transducer probe was centered under the plate- 


footing at a depth of 0.276 inches (7mm). 


To insure proper measurements, the active transducer 
port, PVC tubing, brass coupling and sensing probe were 
thoroughly cleaned with a strong detergent to remove any 
dirt, grease or oil that may have accumulated. After 
rinsing, the active pressure port and extensions were 
carefully flushed and filled with de-aired water. This 
‘procedure was followed to eliminate the possibility of 


entrapped air within the pressure sensing system. 


A static load test was utilized to calibrate the 
pressure transducer. Readings were taken with the 
transducer submerged in still water at a predetermined 
depth. At a second predetermined depth readings were taken 
again. Since the pressure relationship is a linear, an 


average calibration factor of 0.043 psi/mv was obtained. 
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The 5 volt D.C. excitation voltage was supplied by a 
HP17403A pre-carrier amplifier. The output was recorded by 


a HP7402A dual channel strip chart recorder. 


To record the near bed shear velocities at the plate 
edge, a Thermo-Systems Inc (TSI), hot film anemometer system 
was employed. The voltage output of the anemometer is 
proportional to heat loss across the cylindrical filament 
(film) as a result of the fluid flow across the filament. 
The actual velocities were determined from calibration 


curves developed prior to each experiment. 


The calibration was performed by attaching the hot film 
probe to a motorized trolley carriage which rode on a rail 
system above a 120 foot wave flume. The carriage was 


Operated at several different speeds over a predetermined 


distance. The voltage output from the TSI anemometer 
| 


monitor was displayed on a Beckman 200 digital voltmeter and 
recorded. A calibration curve of millivolts versus ft/sec 
es thus developed. Temperature corrections based on 
Reynolds Number similarities were employed to compensate for 


the difference in water temperatures between the wave flume 


,and the test tank. 
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The upward plate velocities (VPU) and the downward 
plate velocities (VPD) were determined from the hot-film 


/anemometer velocity record. 


The scour depth readings were recorded manually by 
observing a marker attached to a stanchion affixed to the 
plate traveling across a stationary rule. Horizontal 


‘deflections were minimized by extending a guide tube over 


l 
| 


the square tubing shaft welded to the plate surface. 


Twelve experiments were conducted at various 


‘combinations of period (T), amplitude of oscillation (A), 


iT) ). 


and plate dimensions (area (AP); weight (W); and thickness 
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CHAPTER IV 
EXPERIMENTAL RESULTS 


A total of 14 experiments were conducted. Two 
preliminary experiments were performed to observe the 
effects of the transducer probe position under the plate. 
As one might expect, the highest negative pore-water 
pressures were recorded at the center of the plate. 
Therefore, the center of the plate was chosen for the probe 
placement. It was the goal of these experiments to obtain 
a general understanding of the mechanisms involved in the 
phenomenon of pumping-erosion. The experiments were also 
designed to investigate whether negative pore-water 
pressures and possible fluidization of the foundation bed 
played a significant role in pumping-erosion. Therefore, it 
was necessary to discern the effects of the experimental 
parameters in pumping-erosion subsidence. Following is a 


table listing the controlled variables for each test run: 
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TABLE 1 - Controlled Parameters 


EXP PERIOD AMPLITUDE PLATE PLATE PLATE 

mo. 6 (T) (A) AREA WEIGHT THICKNESS 
Sec - inches he ae 

ae 7) ee sy hae} ica Se Bo 

2 20 1 1 25 75 

3 5 1 1 25 15 

4 20 1 1 25 75 

5 10 1 1 25 75 

6 5 1 i 25 75 

7 20 2 . 80 36 132 

8 10 2 . 80 36 132 

9 5 2 . 80 36 132 

10 20 2 . 80 36 132 

11 10 2 . 80 36 32 

12 5 2 . 80 36 132 


The soil used was an Ottawa sand with a mean grain 
diameter of .22mm. The average bed density of the 
foundation bed (the area directly under the plate) was 3.30 
slugs/ft3. In addition the controlled parameters listed in 
Table 1, the following variables were also recorded: number 
of cycles(CYCS), positive pore water pressures (PPWP), 


negative pore water pressures (PPWS), plate exit velocities 
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(VE), plate upward velocities (VPU), and plate downward 


velocities (VPD). 


Figures 1 through 12 in Appendix A plot the subsidence 
versus time and cycles. Figures 13 through 24 plot the 


subsidence versus the log of time and cycles. 


Figures 25, 26 and 27 are group plots of subsidence 
versus cycles and Figures 28, 29 and 30 are group plots of 
subsidence versus time for the periods of 5, 10 and 20 
seconds respectively. Figure 31, plots the subsidence 
curves versus cycles (Figure 32 versus time) for the first 
six experiments (1 sf plate). Figure 33, plots the 
subsidence versus cycles (Figure 34 versus time) for the 


last six experiments (2 sf plate). 


Figures 5A through 5L show typical strip chart output 
records for the 12 experiments. The pore water pressure 
record output is recorded on channel 1 and the hot film 


anemometer output record on channel 2. 


As noted in the experimental record, in Appendix B, the 


output of the hot film anemometer became unreliable as the 


| 


water temperature increased above 860F. Thus VE was not 


considered in the analysis described in Chapter V. 
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CHAPTER V 


DISCUSSION OF RESULTS 


As one reviews the results depicted by the figures 
(Appendix A), certain trends appear. The first observation 
is that the experiments that produce the most rapid 
subsidence have the larger plate areas, greater weights, 
larger amplitudes of oscillation, and shorter periods. 
Secondly, the experiments that subside more rapidly also had 


higher negative pore water pressures and thus higher 


liquefaction indexes as derived by Sleath (12), where the 


point of fluidization is regarded as : 


gE fgaste'h 


og oy 0 


water density 


te) 
iT} 


po = bed density 


6p = pore-water pressure through the transducer depth dy 
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Failure of the bed could also result from too large a 


horizontal gradient : 


1 Pins ae “20 tan 6 
Pgh. p 


where 9 = internal angle of friction 


“As the plate is uplifted, as depicted in Figure 6, 
negative pore water pressures are developed locally in the 
foundation soil directly under the plate. The data indicate 
that the bed is indeed fluidized. Assuming that the pore- 
water pressure is incompressible and that the grain skeleton 
is compressible and deformable, the action of local 
fluidization of the bed would than cause distortion of the 


upper grain layers (see Figure 7). 


Some of the loose grains were observed to be carried 
into suspension by the turbulent flow of fluid rushing in 
under the plate in response to the local negative pressures. 
The strip chart record reflects the velocity response to the 
pressure fluctuations which are in agreement with the 


Bernoulli equation : 


p=-}pV2- gx+C 
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During the downward cycle of the plate (Figure 8), the 
loosened grains as well as the sediment in suspension (as a 
result of the uplifting of the plate) can then be easily 
transported (Figure 8) at a much lower shear velocities (u*) 
then would be required for normal initiation as predicted by 
the Shields’ diagram (6). An interesting observation, 
derived from the pressure and velocity records, is that the 
slower periods have more pressure fluctuations and thus more 
velocity fluctuations but with lower negative pore water 
pressures than do the shorter period cycles which have 
greater negative pore water pressures but with a faster 


Change in pressure and velocity with fewer fluctuations. 


By careful observation three interesting details were 
discovered. First, a small ridge of sand developed along 
and under the perimeter edge of the plate. The second 
observation is that a trough was developed immediately 


adjacent to the perimeter edge of the plate (Figure 9). 


Thirdly, slope instability and sloughing was observed during © 


the upstroke of the plate. These seemingly minor details 


may play a significant role in the pumping erosion process. 


The small ridge under the perimeter of the plate 
probably results from the outflow of water and transported 
material attempting to exit from the underside of the plate 
during the last fraction of a second before plate impact. 


This small ridge of sand may account for the large negative 
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pore water spikes at the end of each cycle (Figure 10) in 
one of two ways. First, if the foundation bed is considered 
to be elastic, the negative pore water pressure may be due 
to plate-soil bounce with the small ridge acting to seal the 
plate so that the pressure response is localized. The 
second possibility, is that as the plate impacts on the 
foundation soil, the plate undergoes slight flexure and 
causes the negative pressure excursion. Again, the sand 
ridge would tend to seal off the plate to cause the event to 


be a local effect. 


The trough is caused by the jetting of water from under 
the plate (Figure 9). During the fee of the plate, 
the lower portion of the slope became unstable and sloughed 
onto the aforementioned trough. The entire slope then 
cascaded down to form a new equilibrium slope. The slope 
failure was probably due to the pressure/velocity field near 
the edge of the plate. As the water rushes in, responding 
to the negative pressures under the plate, sand grains at 
the base of the slope could be enveloped by the turbulent 
flow causing a slope failure. An alternate explanation, is 
that there is flow out of the slope in response to the 
pressure fluctuation, which would also cause the slope 


failure. 
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In any event, the sand from the slope failure, which 
should be large grains due to preferential transportation 
Brinnowing), loosely fall into the trough only to be 
expelled and transported during the downward stroke of the 


im plate. 


Topographical plots of the sand contours for the 12 


experiments are in Appendix C. Of particular interest is 
the dynamic angle of repose of the scour mound slopes as 
listed in Table 2. 

TABLE 2 DYNAMIC ANGLES OF REPOSE 


EXP. NO. DYNAMIC ANGLE OF REPOSE 

1 21 

2 24.3 
3 24.9 
4 24.2 
5 23.4 
6 29.38 
7 2353 
8 Zeieie 
9 231.8 
10 (BY 
lig 22.4 


WZ eeT 


iw ,soll Ss? agete cute aoe) Bie Bnd” Jineve % 
inogeuent berénew tap oF sup antag 
od vil fyye \? ay cat. ity yfaaodl } 


> Niswhwod et PPPh Gereaganew 


® 


votneo baa? se %ereselo lentes 
-lujft tao 92 SD whbeeeok at) SRE 


a4 nuom ‘iD? si} 76 gedaet Yo BIpRRG 


WUS3h 19 250A SNARE 
hm — 
Se0ad MH) FA BIMANYG 


——_ 


Figure 11, ranks the experiments by the average initial 
liquefaction index. The ranking is in accordance with the 
observed subsidence rates (i.e. the tests with more rapid 
subsidence have higher liquefaction indexes). A 
relationship between the period, amplitude, plate weight and 
plate area ( A2Ap/WT2 ) was then developed and ranked 
(Figure 12). Curiously enough, the ranking was the same as 
the liquefaction index except for the positions of 


experiments 2 and 4, which were reversed. 


Although some general trends were developed from the 
plotted data and figures, a more sophisticated numerical 
analysis is required to ascertain the relationships between 


the experimental parameters. 


MULTIVARIATE PROCEDURES 


To determine the relationship and relative importance 
of the variables the Statistical Analysis System, or SAS, 
was employed. The SAS techniques were used to analyze the 
data because of the statistical integrity and ease of use. 
All of the statistical methods employed are explained in the 
SAS User's Manual : Statistics (ll). The data from all 12 
experiments were compiled into 234 observations, using 15 
variables. The dependent variable was subsidence (SUBSID). 
The 14 independent variables were : period (T), amplitude of 
oscillation (A), plate area (AP), plate weight (W), plate 


thickness (TH), bed density (BRHO), kinematic viscosity 
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(VISC), number of cycles (CYCS), negative pore water 
pressures (PPWS), liquefaction index (LI), positive pore- 
water pressures (PPWP), plate upward velocity (VPU), plate 
downward velocity (VPD) and plate exit velocity (VE). Table 
3 summarizes the range of values for each of the variables. 
A list of all the variable values can be found in Appendix 


B. 


The SAS procedure, PROC CORR, was first used to 
determine simple linear correlations between the variables 
in the entire data record and are listed in Table 4. Due to 
the high correlation between the area of the plate (AP) and 
the plate weight (W) and plate thickness (TH), W and TH were 
deleted form further analysis. Also LI was deleted from 
analysis due to the high cross correlation with PPWS. 
Furthermore, PPWP was not considered for further analysis 
because of its low correlation (-0.016) with subsidence. 
This is an interesting event as data suggests that positive 


pore-water pressures have little to do with subsidence. 


The next step was to determine the strength of the 
linear relationship determined by computing R2 , the 
square of the multiple correlation coefficient (R-SQUARE in 
SAS). R2 also known as the coefficient of determination, is 
interpreted as the portion of the variability that has been 
accounted for by the regression formula. The closer R2 was 


to 1, the better the model equation fit the data. However, 
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a large R2 did not necessarily mean that the proper model 
was chosen for the data. Therefore, R2 should be used with 
caution, since it is always possible to make R2 large by the 
addition of more variables to the model. The R-SQUARE 
routine was configured to give R2 and C(p) values for all 
possible combinations of 3, 4 and 5 variables. The results 
are listed in Appendix D. C(p), Mallow's statistic, is a 
criteron related to the mean square of a fitted value (7). 
As C(p) attempts to indicates the bias of the prediction, a 
small C(p) value is generally sought. The "best fit" model 
is determined by when C(p) first approaches p on a plot of 


C(p) versus p. 


Next, the remaining 10 variables were incorporated into 
a SAS routine called STEPWISE. STEPWISE is a model building 
program which chooses variables that are to be included into 
the model. Ideally, a variable selection process should 
provide for any given number of independent variables, a 
subset whose equation of predicted values produce the 
minimal residual sum of the squares. In PROC STEPWISE 
variables are added or deleted one at a time, depending on 
the Forward or Backward option selection, until all 
coefficients remaining in the model are statistically 


Siemiticant (P < 0.0001). 
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The STEPWISE procedure, using the backward elimination 
method, gave results that are suspect. It indicated that 
the peeranodel was produced when all the variables were 
significant at the P < 0.001 level when the variable AP 
(plate area) was removed. This result is contradictory to 
the forward selection process which adds variables in the 
order of there significance. The forward selection process 
ranked the variables in order of significance in the 


following manner: 


TABLE 5 FORWARD SELECTION RANKING 


Variables C F PROB>F R2 

GXCS 185.049 Sele. 754 tHa0 4000) 0.6220 
VPD 28.146 143.304 0.0001 0.7667 
VISC 20.663 82842 | 020033 0.7754 
AP 16.634 S737 OL0174 0.7808 
T 14,272 4.209 0.0414 0.7848 
PPWS 11.090 SOSIt F 5c0R0250 0.7895 
A 97307 oe 1G tm O05 3% 0.7930 
VPU 9.000 Mes. 72050302 0.7951 
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The significance of the ranking is apparent in the 
manner in which the forward selection was constructed, in 
that it added the variables in the following order : CYCS, 
VPD, VISC, AP, T, PPWS, A and VPU. Each additional variable 


improved the overall fit of the model. 


The strength of a linear relationship, an indication of 
how well the model fits the data, is indicated by the R2 
value. The F-value is used by SAS to test the null 
hypothesis, which tests the linear association between Y and 
Meee. ... XK . If the F-value is not significant (PROB F 
= 0.0001), for a system on "n" degrees of freedom (DF), the 
null hypothesis is accepted indicating that the model does 
not explain a significant portion of the variability. That 
is to say, as the F-value increases a greater portion of the 


variability is explained by the model (i.e. the spread of 


the model is lessened). 


The PROB >F indicates the significance of the model 
(or the strength of the variables as in the R2 output). The 
lower the value of PROB>F , the more significant the model . 


A p-value of 0.0001 is highly significant. 


The strength of the individual variables in indicated 
by, T and PROB>T, which denote the same inference as F and 


PROB>F. 
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Based on the above results, PROC REG was performed on 
the 8 remaining independent variables. Variables A 
Panel itude) and VPU (plate upward velocity) were then 
removed and PROC REG was run again. The results are shown 


in Tables 6 and 7. 


Two variables, Fl and F2 were then created from the 


combination of the 8 variables where : 


Fl 


{(144 PPWS VPU AP/( g2T))}1/3 


Re 


{(1/2VPD2 A W )/(144 VISC2 g2 T2 )}1/2 


F3 and F4 were created from the log of Fl and F2 taking 
into account the number of cycles. PROC CORR, PROC RSQUARE 
and PROC REG were then utilized on Fl, F2 F3, and F4. The 


results are listed in Tables 8, 9 and 10. 


Plots of the models (Figures 13, 14 & 15) listed in 
Tables 6, 7 and 10, graph the predicted values from the 


models against the actual values. 
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Multivariate Results 


The SAS multiple regression was used to determine if 
the dependent variable Y (subsidence) was linearly related 
to two or more of the independent variables Xl, X2,...Xk . 


The solution takes the form of : 
Y= a+ DIX] + b2X2 +...+ bkXk 


where "a" denotes the intercept term when X's equal 
zero. bl, b2,...bk are the partial regression coefficients 


Onexl, X250..Xk , respectively. 


Below is listed a summary of MOD1, the 8 variable 
model; MOD2, the 6 variable model; and MOD3, the 2 variable 


model containing Fl and F4. 


TABLE 9 MODEL SUMMARY 


SUM OF . 
MOD # VARS R2 F-VALUE SQ. RESIDUALS 
1 8 OR 795r 109 1.78 
2 6 0.7895 142 12.109 
3 Zz 0.6618 226 19.43 


* 2 terms containing 8 variables 
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From Table 11, the 6 variable model was selected as the 


best model where: 


fea cu0er -) 0 11987 > 0.7103 AP + 000536 CYCS + 
LS ol 97 PPWS. + 21645259. VISC + 0.16325 VPD 


The determination of the "best" linear model is not as 
important to this research as is the significance of the 
variables. From Tables 5 and 7, it can be seen that the two 
most significant variables are the number of cycles (CYCS) 
(F = 381.754, T = 18.39, p-value = 0.0001) and downward 
Dilace velocity (VPD) (F = 143.304, T = 10.225, p-value 
=0.0001). The remaining variables listed in their order of 
significance are : viscosity (VISC), plate area (AP), 


Negative pore water pressure (PPWS) and period (T). 
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CHAPTER VI 


SUMMARY AND CONCLUSIONS 


The main objective of this research was to determine 
which mechanisms were involved in the phenomenon of pumping- 
erosion. The secondary objectives, involved the evaluation 
of the relative importance of each of the 14 independent 
variables, and in particular, the significance of negative 


pore-water pressures, 


Of the original 14 independent variables, 6 variables 
(cycles (CYCS), downward plate velocity (VPD), viscosity 
(VISC), plate area (AP), negative pore water pressure 
(PPWS), and period (T)) were considered significant in the 
model constructed utilizing SAS. The number of cycles 
(CYCS) and the downward velocity of the plate (VPD) were 
highly. significant and is in agreement with the accepted 


understanding of pumping-erosion and scour. 


The result of negative pore water pressure (PPWS) 
remaining significant throughout the analysis supports the 
hypothesis that negative pore water pressures are a 
contributing factor in pumping-erosion. This is 


understandable if the experimental results of the 
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Liquefaction Index (fluidization index) associated with the 


negative pore water pressures is accepted. 


This leads to the conclusion that the fluidization of 
the foundation soil, due to the negative pore water 
pressures developed during the upward cycle of the plate- 
footing is a precursor to the scour developed during the 
downward cycle of of the plate-footing and significantly 


contributes to the overall pumping-erosion subsidence . 


As with many initial investigations, this research 


poses more questions than it answers. Such as : 


1) If negative pore-water pressures and the 
subsequent fluidization are significant to the 
phenomenon of pumping-erosion, then what is the 
distribution of the negative pore water 


pressures under the footing ? 


2) What is the distribution of velocities under 


the plate footing ? 


3) Can the footing be designed to alter the 
pressure and flow nets under the footing to 


reduce subsidence ? 
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4) What effect does the shape of the scour mound 


have on the negative pore water pressures ? 


5) How would the results vary in other 


cohesionless and cohesive soils ? 


6) What are the effects on pivotal oscillations 


as opposed to vertical oscillation ? 


RECOMMENDATIONS 


To obtain a better understanding of pumping-erosion, 
more research is needed in this area. If the results are 
valid, the effect of negative pore water pressures may 
affect all dynamic structure-marine soil interactions 


including submarine pipelines. 


For purposes of further research the following 


recommendations are proposed : 


1) Employ hydraulic or pneumatic rams to load the 
| soil at a constant rate with a constant 


amplitude. 
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2) Develop a method to record the dynamic 


formation of the scour mound slope. 


3) Conduct the experiments with different grain 


sizes of cohesionless soils. 


4) Conduct the experiments with cohesive soils. 
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Alvin Eugene Grimmig Jr. was born on May 24, 1951 in 
“Miami Beach, Florida of Alvin and Nora Grimmig Sr. He 
married Miss Debra Kay Rademacher on May 6, 1972. He 
enlisted in the Navy in June 1972 and received a commission 
in Civil Engineer Corps in 1978. He received a 
baccalaureate degree in ocean engineering from the 
University of Washington located in Seattle, Washington in 


December of 1978. 


Permanent address: 


LT Alvin E. Grimmig Jr. 
c/o Barbra Rademacher 

555 Filmore Ave. Apt.101 
Cape Canaveral, FL 32950 
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Subsigence vs Time 
EXPERIMENT 4 T= 20 secs A= 0.8 inch Ap= 1sf Wp= 25ibs 
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Subsidence vs Cycles 
EXPERIMENT 4 T= 20 secs A= 0.8 inch Ap= 1sf Wp= 25lbs 
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FIGURE 4 
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Subsidence vs Time 
EXPERIMENT § T= 10 secs A= O.8 inch Ap= Isf 
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Subsidence vs Cycles 
EXPERIMENT 5 T= 10 secs A= 0.8 inch Ap= Isf 
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FIGURE 5 
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Subsidence vs Time 
EXPERIMENT 6 T= 5 secs A= 0.8 inch Ap= isf Wp= 25lbs 
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Subsidence vs Cycle 
EXPERIMENT 6 T= 5 secs A= 0.8 inch Ap= 1sf Wp= 25ibs 
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FIGURE 6 
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Subsidence vs Time 


EXPERIMENT 7 T= 20 secs A= 1.0 inch Ap= 2sf Wp= 36ibs 
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Subsidence vs Time 
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